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Abstract: We study the scattering of a massless scalar field in a generic Kerr back¬ 
ground. Using a particular gauge choice based on the current conservation of the radial 
equation, we give a generic formula for the scattering coefficient in terms of the com¬ 
posite monodromy parameter a between the inner and the outer horizons. Using the 
isomonodromy flow, we calculate a exactly in terms of the Painleve V r-function. We 
also show that the eigenvalue problem for the angular equation (spheroidal harmonics) 
can be calculated using the same techniques. We use recent developments relating the 
Painleve V r-function to Liouvillc irregular conformal blocks to claim that this scat¬ 
tering problem is solved in the combinatorial sense, with known expressions for the 
r-function near the critical points. 
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1 Introduction 

Black hole scattering theory is an important problem from the general relativity and 
string theory perspective. It serves as a way of testing linear stability [1], superradiance 
[2], black hole entropy [3] and relaxation times in AdS-CFT [4], For the rich story be¬ 
hind it see [5, 6]. Scattering is also important to study astrophysical problems related 
to the detection of gravitational waves and a possible mechanism for short gamma ray 
bursts, for example. See the introduction of [7], references therein and [8, 9] for discus¬ 
sions. It is also a daunting problem from the mathematical side, involving unknown 
properties of new special functions [10]. Most of the previous study in four dimensions 
on the subject relied on numerical analysis, either by direct numerical integration or 
by expansion of the solution in terms of known special functions [11-13]. Also, decou¬ 
pling limits like the near-extremal cases were taken directly on the metric (see [14] for 
applications in Kerr-CFT). 

In this letter we build on previous analytical work, a reverse Hilbert-Riemann prob¬ 
lem [15], based on isomonodromy equations [16] to relate the scattering coefficients for 
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a massless scalar field in four-dimensional Kerr background in terms of the Painleve V 
transcendent r-function [17] and show that the problem can be solved in the combi¬ 
natorial sense using the expasions of the r-function near the critical points. We also 
show that the angular equation can be solved implicitly using the same techniques, in 
a new way of extracting eigenvalues for the spheroidal wave functions [18, 19]. The 
results shown here can be understood as the natural conclusion to the application of 
the monodromy problem to solve scattering problems as applied to black holes [20-23]. 
The result follows from a careful consideration of the “linearization” of the Painleve V 
equation in the sense given by Miwa, Jimbo et al. in [24-26], and of the time-reversal 
symmetry. These two points allows us to place boundary conditions such that the two 
linearly independent solutions of the radial equation correspond to the two solutions in¬ 
volved in the connection problem of the Painleve V equation [27]. Recent developments 
on the relation between Painleve transcendents and the AGT conjecture [28, 29] show 
the striking connection between black holes and black hole scattering to Conformal 
Field Theory [30]. 

This letter is organized as follows. In Section 2 we treat the radial equation for the 
massless Klein-Gordon field in a generic four-dimensional Kerr background. We write 
the equation in a canonical confluent Heun type and discuss the time-reversal symmetry. 
In Section 3 we discuss how the scattering problem for the radial equation is connected 
to the connection problem of the confluent Heun equation. We describe the (extended) 
monodromy data of the confluent Heun equation, including the Stokes parameters, and 
use the time-reversal symmetry to parametrize the scattering coefficient in terms of the 
monodromy parameters. In Section 4 we discuss how to extract these parameters from 
the Painleve V r-function, using the isomonodromy flow. In Section 5, we outline a 
procedure to extract the eigenvalues for the angular equation using the same structure. 
Finally, in Section 6 we list the results of [17] which give a combinatorial expansion for 
the Painleve V r-function near to = 0, which would correspond to our extremal limit. 
We close with a short discussion and future directions in Section 7. 

2 Klein-Gordon in Kerr 

It is known that the massless Klein-Gordon equation in a Kerr background is separable 
[31], i. e., the solution can be written as = e~ lbJt e imcl> R{r)S{6). The angular 

part S{6) is solved in terms of oblate scalar spheroidal harmonics, whose eigenvalues 
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(the separation constant C'y m ) can be approximated for different asymptotic regimes, 
including low and high frequency limits [18, 32], and will be the object of study in 
Section 5. We will focus for now on the radial part 


d r {Q{r)d r R{r)) + 



R(r) = 0 , 


where Q(r) and W 2 are given by 


Q(r) = r 2 — 2 Mr + a 2 = (r — r + )(r — r_), 
W r = u(r 2 + a 2 ) — am. 


( 2 . 1 ) 


( 2 . 2 ) 


The form which will be useful for the analysis has 


0 = 2 iu(r — r_), and y(z) = (r — r_) 9 °^ 2 (r — r + ) 9t °^ 2 R(r), (2.3) 


with the parameter 

t 0 = 2 iu)(r + — r_). 

The radial equation is written in this parametrization as 

+ pWt + q(z)y(z) = 0, 

az z az 


P(z) = 


1 - 9 to 1 - 9 0 , . 1 , c 0 c to 


+ 


9(*) = ~ + - + 


z — to z ' "' 7 4 0 z — to' 

which is known as the confluent Heun equation, with parameters given by 


9o = 2 i 


,uj(r_ + a 2 ) — am i oo — Vtm 


r_ — r + 


2vr T_ 


0 to = 2i 


,co(ri_ + a 2 ) — am i u — ttm 


r _ r _ 


2vr T, 


t 0 C 0 — Cl,m + — ((1 — to)6o + 0t o ), toCt 0 — —C^ m — -(00 + (1 + to)9t 0 ), 


(2.4) 


(2.5) 


( 2 . 6 ) 
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where we defined the temperatures T± and angular velocities Q± for the outer and inner 
horizons. For future reference, we also define 

$oo = 1 + 2(co + c to ) = 1 — AiMu = 1 — 9 0 — 9 to . (2-7) 

Also, eq. (2.1) has real coefficients, so if R{r ) is a solution, so is R{r)*. We then define 
the radiation flux as 


j = -iQ(r)[R(r)*d r R(r) - R(r)d r R(r)*}. (2.8) 

One should bear in mind that real r means z* = —z, because of (2.3). 

3 The Fuchsian System 

Writing the solutions for (2.5) in terms of Frobenius expansions at z — 0 or z — to 
is straightforward. On the other hand, extracting the scattering coefficients from the 
solutions is not straightforward in the sense that not enough information is known about 
the function defined by the ODE. To wit, we are interested in the so-called connection 
problem: how a solution of (2.5) with known behavior at one singular point, say yt 0 (z ) 
at z = to, is related to the pair of solutions with known behavior at another singular 
point, say, y^(z) at z = oo. So, for example, in 

Vto(z) = Cy+(z) + Dy~(z), (3.1) 

C and D are the connection coefficients. If yt 0 (z) represents a purely ingoing wave 
at the outer horizon z = to and y^(z) ingoing and outgoing waves at infinity, then 
the connection coefficients are indeed the scattering coefficients. As we will see below, 
solving the connection problem allows us to solve the scattering problem. 

In [16], the authors proposed a method based on the Schlesinger equations to 
extract the scattering coefficients. It relies on viewing the ODE (2.5) as the equation 
satisfied by the first row of the 2x2 matricial Fuchsian system 

= A(z)$(z), (3.2) 

and to extract the connection coefficients from the monodromy data of the Fuchsian 
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system. The monodromy data for the system at hand is defined as follows: suppose 
A(z) is meromorphic and let {^} be the isolated singular points of the system (4.1). 
The point at z = oo is studied by the substitution u = 1 /z. The solution <3>(z) will 
also cease to be analytic at these points and behave singularly there. Since <f>(z) = 
®(e 2m (z — Zi) + Zi ) is also a solution of (4.1), then we must have <E>(z) = <P(z)Mj, with 
Mi a constant matrix associated with a circuit around Zi - the monodromy. 

The form of Mi depends on the behavior of A(z) near the singular point Z\. If A(z) 
diverges like (z — ^) _1 , then the singular point is called regular and Mi is conjugated 
to diag(e 27ria * , e 2?ua * ), where af are the solutions of the indicia] equation near z, L . This 
form of Mi stems from the leading behavior (z — Zi ) a * of the solution near z t , a branch 
point for <f>(z). We will for now on assume that the differences af —a are not integers 
for the radial equation. If, however, A(z) diverges like {z — Zi)~ 1 ~ ri , where the integer 
rj > 0 is the Poincare rank, the solution will have an essential singularity and will ex¬ 
perience the Stokes phenomenon-, the leading term in the Frobenius expansion displays 
an exponential behavior which diverges or converges depending on which direction the 
limit z —» Zi is taken in the complex plane. In order to describe it, we will specialize 
to our case. In equation (2.5), the matricial system has A(z) with the structure 

-4W = ^3+ — + (3.3) 

2 z z — t 

where there are two regular singular points, at z = 0 and z — t, and an irregular point 
at z = oo with Poincare rank 1. We also have the usual Pauli matrix cr 3 . The elements 
of Ao and A t in (3.3) are independent of z and are obtained from the ODE (2.5). This 
identification has extra freedom because the Fuchsian system has more free parameters 
than the ODE. We thus pick the gauge choice 

Tr Ai = 9i, det A; = 0, Tr a 3 (A 0 + A t ) = -9^, (3.4) 

which, along with the ODE parameters (2.6), fixes the elements of Ai up to conjugation 
by e urT3 . The epithet gauge choice will become clear in the next Section. 

In order to describe the system at the irregular singular point z = oo, one defines 
the sectors [27, 33] 

Sj = {z e C| (2 j - 5)| <&vgz < (2 j - 1)^|, (3.5) 
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for j G Z. On each Sj, one defines the asymptotic behavior for 'h(^) = z 2^ e ° +et ' > ^(^z) 
based on the gauge choice (3.4) for A(z), 


T^)^. = (1 + O^ 1 )) exp(±zcr 3 )z 2 0OO,T3 


(3.6) 


These canonical solutions are connected by the Stokes matrices T J+ i(z) = T j(z)Sj 
and, given that ^j +2 (e 2m z) = ^ j (z)e~ md ° oCr3 , we have S j+ 2 = a, and 

the Stokes matrices are all determined by any consecutive 2, which have the following 
structure 


S V = 



S2j+1 ~ 



(3.7) 


and the S 2 j and s 2 j + i are called Stokes parameters. Local solutions of the Fuchsian 
system at z = oo are defined in the Riemann surface of the logarithm, and these 
parameters describe how the solution (3.6) changes as one changes the leaf of the 
surface. The (extended) monodromy data for (2.5) is completely determined by the 
monodromy matrices M 0 and M t , at the regular singular points 0 and t, respectively, 
and the Stokes matrices So and Si, which determine the (formal) monodromy matrix 
at oo in sector Sj by 

= SjSj^e™ 6 ™" 3 . (3.8) 


Knowledge of the and of two consecutive Stokes parameters, say, Si and s 2 is sufficient 
data for solving the connection problem. We define the connection matrices E t by 


Mi 





i = 0,1, 


(3.9) 


and those are given in terms of 6^, a, s i and s 2 in [27, 33]. Note that, because 
MooM to Mo = 1, we have that the combined monodromy parameter 2 cosncr = Tr M 0 M t 
satisihes 

2 cos 7rcr = 2COS7T0OO + e” r(9oo s 1 S 2 , (3.10) 

and that M. ^ = E^e inaa3 E^. The connection matrices Ei can be seen - when properly 
normalized by det £7=1 - as the matrices that implement the change to the “natural 
basis” of solutions of the ODE at Zi where the monodromy Mi is diagonal. Explicit 
formulas for E- t in terms of monodromy data are given in [33]. We use the remaining 
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gauge symmetry A(z) —> e urT3 A(z)e~ urT3 so that both entries in the first row of <F(z) have 
equal but opposite current density normalization (2.8), defined as j ± . So we have the 
entries oc i?+ and [$(z)]i 2 = oc being the normalized incoming 

and outgoing waves near infinity. Because the are conserved, this choice also means 
that the entries will be normalized near the outer horizon to, up to a phase. Therefore, 
the entries of the connection matrix 


Ec »to = E^E ^ 1 


(3.11) 


will give the explicit scattering coefficients, allowing us to write T and 1Z such that 


,+ 


n 

— 


u t = —u^ , 


(3.12) 


and by construction the time-reversed wave is equal, up to a phase, to the entry [<E>(z)] 12 
near to- We notice that Det_E' oo t 0 = 1 implies \1Z\ 2 = 1 — |T| 2 , which is the condition of 
flux conservation. Now, by substitution of the diagonalized forms of the monodromy 
matrices in Tr = Tr Mq 1 , and using (3.12), we obtain the formula 


m 2 


sin 7rcr sin TT0 to 

sin f (0 O + a - 0, a ) sin f (0 O - a + 0 h ) 


(3.13) 


which gives the coefficient in terms of the composite monodromy parameter a. In order 
to compute it, we will need the structure of isomonodromy deformations, outlined in 
the next section. After embedding this system in the isomonodromy formulation, we 
will be able to give explicit coordinate dependence of the wave functions. 


4 The Isomonodromy Method 

The idea of considering the Fuchsian system (4.1) instead of the scalar ODE came in the 
context of the Riemann-Hilbert problem, since a simple counting argument [35] shows 
that the generic scalar ODE does not have enough parameters to realize all possible 
monodromy matrices. The Fuchsian system, however, has too many parameters, so 
there is a family of A(z) in (4.1) with the same monodromy data. The description of this 
family can be understood physically: suppose we are given an invertible solution < h(^) 
of (4.1) with some singular points. Then we can understand A(z) = [■Jr < T , (^)]4>(^) —1 as a 
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“gauge potential” describing non-abelian charges, dipoles or multipoles at each singular 
point. The physical configuration described by the gauge potential is invariant under 
generic gauge transformations A(z) —» U(z)A(z)U(z)~ 1 + [-jpU(z)]U(z) -1 , meaning 
that the monodromy data (non-abelian holonomies) are also invariant under these 
transformations. 

In the case at hand (2.5), with A(z) of the form (3.3), we supply a t-component to 
this gauge potential, 

B(z) = (4,1) 

z — t 

in such a way that the “field strength” dtA — d z B + [A, B] is zero if the Schlesinger 
equations hold 

8A 0 1 r . , , dA t 1 , . . , 1 r , , . , . 

~df = t [At ’ Ao] ’ ~dt = — t ^ A ^ ~ 2 Atj 0-3 ’ (4 ' 2) 

The monodromy data of the Fuchsian system (4.1) will be independent of t if Ao(t) 
and At(t) in (3.3) satisfy the equations above (4.2). For the relationship between the 
Schlesinger equations and the theory of flat holomorphic connections we refer to [24- 
26, 35], 

The system (4.1) has some direct conserved charges. The gauge choice (3.4) is 
justified now because these are the constant of motions under Schlesinger evolution. 
With these constraints, the equation satisfied by any element of the first row of $(z) 
in (4.1), with A(z) given by (3.3), is of the form 


p(z) 


d 2 y 

dz 2 


+ P(z ) 


dy_ 

dz 


+ q(z)y = 0, 


1-flo 1-0* 1 

z z — t z — A’ 


q(z) 


1 c 0 

-7 + — + 

4 z 


C t 

z — t 


+ 


d 

z — A ’ 


(4.3) 


where Co, Ct, A and /i are functions of the entries of A(z). Compared to (2.5), the 
system above has an extra singularity at z = A. However, this is an apparent one - its 
monodromy is trivial, as the indicial equation has exponents 0 and 2 and there is no 
logarithm behavior due to the relation between fi, A and t, 


d 


2 


0n-l 


+ 


0*-l 
A — t 


- 8 


C 0 

,< + T + 


c, 

A — t 


1 

4 


A 


(4.4) 















The system (4.2) is better known in another clothing, as 


/,n _ [Ao(t)]n[A t (t)] 12 _ d 0 + 6t — 600 — (2/i — 1)(A — t) 
[At(t)]ii{A 0 (t)]i2 9q + 9t — doo — (2p — 1)A 


(4.5) 


can be checked to satisfy the Painleve V equation. These are part of the Painleve tran¬ 
scendents family of differential equations: non-linear second order differential equations 
whose solutions do not possess movable branch points [35]. These define new special 
functions, with applications in integrable systems, random matrix theory and confor¬ 
mal field theory - see [17, 36] for references in those applications. I 11 [27], asymptotic 
expressions for the Painleve V system were given in the guise of the T-f u n otion 



0CO 11 dodt (4 - 6) 


in terms of the monodromy data. The r-function is the most natural isomonodromy 
invariant that can be defined, with a clear interpretation as a generating functional in 
quantum field theory applications [37-39]. 

As the problem now stands, knowledge of the r function solves the system com¬ 
pletely. We pick suitable initial conditions for p and A in (4.3), and from those recover 
the A 0 and A t , which then will set the initial conditions for the r function. By choosing 
the initial conditions 


&t — @t 0 ~ 1) A(£ 0 ) — to, fJ>(t 0 ) — -r——- 

dt n ~ 1 


(4.7) 


based on the parameters of the initial Heun equation (2.6) we have 



(4.8) 



which can be formally inverted to yield the Stokes parameters Si and s 2 , and hence a. 













We note that our definition (4.6) differs slightly from the one in [27] by 



(4.9) 


One also notes that (4.8) allows us to interpret the r-function as the generating function 
that implements the canonical transformation between the non-trivial monodromy data 
si, s 2 and the canonically conjugated parameters of the ODE t 0 , c to . 

The structure outlined here gives another interpretation for the r-function. One 
of the most fundamental results from Miwa, Jimbo et al. on the Schlesinger system 
[24-26] - see also [35] - was the Hamiltonian structure of the isomonodromy flow. This 
structure stems from the symplectic structure of flat holomorphic connections (see [16] 
for a description), but can be intuitively understood from the algebraic constraint (4.4). 
One first notes that C 0 and C t are not independent, they are related to the asymptotic 
behavior at z = oo, via 60), 



(4.10) 


Varying t will now tie the variation of A and /j in a symplectic system, in which Co 
can be thought of as the Hamiltonian function. Now, by “solving” the Hamiltonian 
system one means giving a canonical transformation from the coordinates Co, to to 
coordinates where the isomonodronic flow is trivial. But the latter is parametrized 
by the monodromy data, which by construction is invariant under the flow. The r- 
function is then the generating function that implements this canonical transformation, 
as evidenced by (4.8). In [40] a set of Darboux coordinates were constructed for the 
Painleve VI case. It would be interesting to give a parallel of that construction to the 
case at hand. 

4.1 Asymptotic wavefunctions and normalization 

Now, with the embedding of the radial equation (2.1) in the matricial system via (4.7), 
we can work out the specific form of the wave functions in terms of the coordinate r. 
The form is not necessary for our results, but will help comparing with previous work 
[12, 13]. The normalization of the elementary matrix 4*(^) is given by 


det$(z) = z e °(z -to)' 940 " 1 


(4.11) 
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with boundary conditions at z = oo given by (3.6). It is straightforward to verify that 
the Wronskian between the two entries of the first row is given by: 


W(z) = [$'(^)] 11 [$( 0 )] 12 - [$'(^)] 12 [$(z)] n = [A(z)\i 2 det $. (4.12) 


Now, using the form of [A(^)] 12 compatible with (4.3), with conditions (4.7): 


we arrive at: 


W*)]l2 = 


k(z — A) 
z(z — t) 


k 

z 


W{z) = kz e °~ 1 {z - t 0 ) et o~ l . 


(4.13) 


(4.14) 


We remind the reader that the value k is “gauge dependent”, and proceed to normalize 
the asymptotic from the entries of $( 2 ). From (3.6), we have: 


[$oo(z)]n = + 0(z~ 1 )), 

1 1 (4-15) 

[4>oo(^)]i2 = kz2 {9o+et o-i+Ocol-ig-a^i + O^z- 1 )), 

where the limit is taken through the direction argz — (z —>■ +ioo). Following [27] 
we define the Frobenius basis at z — to to be: 


^*0^) - Gt ° ^0 ( z - fo)^" 1 ) ^ + °^ Z ~ ^' 16 ^ 

where Gt 0 is a matrix which diagonalizes A t . Its particular form will not be important 
for us save that for the conditions (4.7) it can be checked to be of the lower triangular 
form. From this we can write the Wronskian-normalized basis for wavefunctions at 

2 = t 0 : 

[$to(^)]n = tf°(l + O(z-t 0 )), [$t 0 (z)]i 2 = ktf 0 ~ 1 {z-t {) ) et ^{l + 0(z-to)) } (4.17) 


Now, we can write the current-normalized wavefunctions in terms of the coordinate 


- 11 - 



r: 


u, 


'c 



As anticipated in the previous section, a suitable choice of k can make the relative 


normalizations equal. The result (3.13) follows. In order to implement the time- 


reversal procedure, one has to check the relative phases and include them explicitly in 
the formula for the connection coefficients. At any rate, the basis at r = r + is different 
from the one considered by [12, 13], which were mainly concerned about generic spin 
Teukolsky equation and did not work the normalization out explicitly near r = r + . 
In the non-zero spin case, one has to solve not only for the monodromy but also for 
the Teukolsky-Starobinski identities [41]. It is an interesting future problem to try and 
extend the results presented here to that case. 

5 Angular Equation 

The general solution for the composite monodromy parameter a using (4.8) (and (3.10)) 
can of course be applied to more generic situations. The associated equation for the 
angular part of the wavefunction S(8) is also of the confluent Heun type. Following 
[19], the scalar spheroidal harmonics are the solutions of the angular equation: 



for which the solution has the following behavior near the singular points x = cos 9 = 


± 1 : 



(1 + x)"*/ 2 , x^-1, 

(l_ x )_ m/2 , 


(5.2) 
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for integer m. It is clear that this behavior will only happen for particular discrete 
values of the constant Cg tTn , which can be fed into (2.1). 

To bring this equation into the canonical confluent Heun form (2.5), we make: 

z = 2auj(x — 1), y(z) = (1 + x) m/2 (l - x)~ m/2 Se, m (x), (5.3) 


and the new parameters are: 


t 0 = -4 auj, 9 0 = -9 to = —m, c 0 = -c U) 


Cl,m + 0 ?UJ 2 

Aau 


(5.4) 


Now, from the discussion above, the existence of a solution of the confluent Heun 
equation with behavior given by (5.2) puts constraints on the connection matrices. In 
this particular case, one of the natural solutions at z = 0 will not “mix” with the natural 
solution at z = 1. This would correspond to the vanishing of one of the constants at 
(3.1). One should not that the second natural solution will of course mix. Also, since 
the parameter 6q and 6t are integers, there will be logarithm behavior for the other 
solution - in the case to —* 0 it would correspond to the associated Legendre functions 
of the second kind. 

It is not hard to see that from this condition the composite monodromy 2 cos na = 
Tr M 0 M t will have to be special: 


a = 2£, £ e Z, 


(5.5) 


as required by the vanishing of the analogue of (3.13) for this case. Now, using (4.8), 
one has a - somewhat formal - solution for the eigenvalues: 


Ce,m = —a 2 0J 2 — m{m — 1) + t — log r(i; m, — m, 2£) 


(5.6) 


t=—4acd 


Since r(t) is also a function of the Stokes parameters si, s 2 , the second equation in (4.8) 
is also necessary to find the Cg. 


6 Asymptotics of Painleve V 

In the spirit of making the analysis self-contained, we copy the relevant formulae about 
the Painleve V r-function from [17]. In the following we use their definition for the 


13 




T-function: 


r(t) = (6.1) 

and assume a and 0* — 6j, a — dj are not integers. The expansion for the tau-function 
is of the form 


f(t, 0) = ^ C({0i}, \cr + n)s n t^ a+n)2 B{{6i}, \a + n; t), (6.2) 

nEZ 

where the irregular conformal block B is given as a power series over the set of Young 
tableaux Y: 

= £ B A ,„({e i },iff)t |A|+l '*i, (6.3) 

A,AtSY 

with coefficients 


Ba„= n 

(*,j)e a 


'4$oo + ^ cr + * 


+ \° + i ~ j) 


•\2 


l^o 2 ) 


+ to - i - 3 + 1 + <r) 


x 


n 

(hj)€V> 


In 

’ 2 ^ oo 


^ + i - i)((|^t - ^ + i - j) 2 - \0 q) 


h 2 M,j)(Xi + fi , i -i-j + l + a) 


(6.4) 


where A denotes a Young tableau, A* is the number of boxes in row i, A' is the number 
of boxes in column j and h\(i,j) — A, + A'- — i — j + 1 is the hook length related to the 
box (i,j) G A. The structure constants C are rational products of Barnes functions 


c«(vu)=n 


G{1 


+ e|cr)G(l + ^6t + ^9o + e\o)G{l + | Ot — + e ^ cr ) 


e=± 


G(1 + ea) 


(6.5) 

where G(z) is defined by the functional equation G{l + z) = T{z)G{z). The parameters 
cr and s in (6.2) are related to the “constants of integration” of the Painleve V equation. 
In our treatment, they are functions of the Stokes parameters Si and S 2 - The parameter 
a is given by (3.10), whereas the expression for s in (6.2) is rather long and involved. 
We will outline the procedure in Section 10 of [33] to compute it. Let M 0 , M t be 
the monodromy matrices for the hypergeometric equations. They are of the form 
Mi = EMe mei<T3 Ei with E t given by equation (10.15) in [33]. They satisfy M t M 0 = 
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e " r ( 6, ° +6ll )diag(e ma , e ma ). We introduce the matrix 


( r(-o-) r(—cr) ^ 

r(i-i( CT -6/cc)) ra-i^+eoo)) 

^in(a— gno)/2 r(cr) T(—a) 

\ r^+floo)) r^a-e^)) J 


( 6 . 6 ) 


whose purpose is to match the diagonal monodromy of the hypergeometric system at 
infinity to the non-diagonal monodromy of the confluent hypergeometric system at 
the origin. It is defined up to a diagonal matrix s _<T3 . We now have the monodromy 
parametrization for E t in (3.9): 


Ei = EiS-^SR-* 3 , (6.7) 

where R 2 = S\/s 2 serves to bring M 0 0 to a symmetric form. With this intricate con¬ 
struction, the parameter s can be read implicitly from M t = M^Mq 1 - or more easily 
from its trace. The parameter s in (6.2) is given by s = — ^(cr + 6 l 00 )s 2 . These expres¬ 
sions are purely algebraic in the sense that they can be computed from the constraints 
M^MtMo = 1 obeyed by the monodromy matrices. The expansion (6.2) is obtained 
from the AGT conjecture [28], relating the Liouvillc conformal blocks to the instanton 
partition function of Af = 2 SUSY quiver theories in four dimensions via a confluence 
limit. 

The asymptotic limit of the Painleve V r-function was also considered in great de¬ 
tail in [27, 33]. For our purposes, in the limit t 0 —>■ 0, the sum (6.2) can be approximated 
by three terms: 

f{t,{e i },a) = e~i e * t &- 2n ^f(t), (6.8) 

/(f) = K( 1 + Bit + Uif 1 -^- 2 ^ + e>(f 2 , t 2{ i±U-2n)))) j (6.9) 

where we assume that 2n < a < 2n + 2. The value of n e Z will be obtained below, 
and we will also assume that g — a — 2n < 1, which can be accomplished by shifting 
n, and the calculation that ensues will be analogous as below. As per (6.2) the leading 


15 








parameters of the expansion of f(t) are: 


k — c({ 0 i}, \o — n), B t = ~(e; - el + e 2 ), (6.10) 

Cl = 8^“i-e)2 ((<, ‘ + e) ' “ S ° 2)r ‘' (6 ' n) 

Now, (4.8) is written as: 

j((9o — 9 to + l) 2 — g 2 ) — toy + 7jto(9 to — 1) = toCt 0 — 9o(9 to — 1) (6-12) 

-toy- + 4j>= 0- (6-13) 

As alluded above, the second equation does not entail information about the initial 
value of the isomonodromy flow beyond the value of to- We will see it as fixing the 
value of s in terms of g. Plugging the expansion for /(f) as above, we find: 


c ‘ = - +(6 - 14) 

Now, back to (4.8) and isolating g, we find: 

Q 2 — (9o + 9 10 ) 2 — 1 + 4 C( m — 2 to + 4 t 0 9 to + 4t 0 —4 -—Bi + 0(t 1+e ), (6.15) 

i - e 

where e = min(p, 1 — g). To first order in to, we then have: 

Q 2 - Go- U) 2 - A9 t0 -—^—^{(9 t0 - l) 2 - 9 2 q + qI) , (6.16) 

£ 2 = -l-16MV + 4C/ m , (6.17) 


where we see that the integer n is chosen to offset the integer part of C/ m . Also, for 
small cu, go is purely imaginary, but the to corrections will introduce a real part, and the 
expansion will have poor convergence. This is verified for the Painleve VI expansions 
using formulas analogue to (6.2) in [17], as convergence of the series improves for finite 
to- Notwithstanding these remarks, we find qualitative agreement with [12, 15] even 
though the small to expansion is not strictly speaking a low-frequency expansion, as 
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other parameters also depend on u. 

At any rate, the expression (4.8) should be valid for generic values of the parame¬ 
ters, and with the results of Section 5 where the separation constant Cy m are also com¬ 
puted using isomonodromy, we can claim the formal analytical solution of the scattering 
problem is presented. We leave for future work a survey of the values for the scattering 
coefficients, as well as expansions for the near extremal and near-sup err adiance 6q ~ 0 
cases. 

7 Discussion 

In this letter, we sucessfully related the scattering coefficients for a massless scalar field 
in a Kerr background in terms of the Painleve V r-function. By choosing the current 
normalization of the ingoing and outgoing waves, we first establish that the expres¬ 
sion (3.13) does give the transmission coefficient directly in terms of the composite 
monodromy parameter a. Also, by relating the problem of finding a to the inverse 
Riemann-Hilbert problem, we are able to give analytic, if implicit solutions for the 
composite monodromy (4.8) and the eigenvalue of the angular equation (5.6) in terms 
of the Painleve V r-function. Also, approximate expressions for small isomonodromy 
parameter - related to low-frequency and/or near extremal case - can be obtained from 
the asymptotics of Painleve V equation in terms of more elementary functions. It is 
our hope that the application of the techniques outlined here and in previous work by 
the authors [16] will help solve many outstanding problems in scattering for black-hole 
backgrounds, like Kerr-(A)dS, and higher spin perturbations. It should be stressed that 
the usefulness of (4.8) goes far beyond the low-energy or the near horizon expansions, 
- both tied to small t 0 - allowing one to study normal modes, and the onset of superra¬ 
diance at 9o ~ 0. One can also use essentially the same arguments given here to study 
scalar perturbations between the inner and outer horizon, providing yet another tool 
for testing for instabilities. We will leave these important directions to future work. 

The result showed here relied on the conservation of current argument, which al¬ 
lowed us to equate the scattering problem to the (mathematically defined) connection 
problem of the radial equation. Although this needs proof, the same argument should 
work for arbitrary spin perturbations, the so-called Teukolsky Master Equation (TME), 
whose radial and angular equations for Kerr background are also separable and yield 
the same confluent Heim type considered here. This should provide a powerful analyt- 
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ical tool to deal with normal modes, generation of gravitational waves by black holes, 
and the linear stability of the inner horizon. 

In a more abstract viewpoint, the results extracted here point to the striking re¬ 
lationship between black hole physics and Conformal Field Theory. Relations for the 
r-function (6.2) were obtained because of recents results concerning Liouville Field 
Theories in the semi-classical limit, where the Painleve V r-function arises as the ex¬ 
pansion of the irregular conformal blocks, which are obtained from a confluence limit 
of (normalized) 4-point functions in Liouville Field Theory (see [17] and [16] for de¬ 
tails). The 4-point function itself arises as the Painleve VI r-function, and in the Black 
Hole case in the Kerr-de Sitter conformally coupled scalar field scattering [42], The 
confluence limit of the Painleve VI is the same as the zero cosmological constant limit. 
This can be seen somehow as an extension of the ideas of [30], where conformal sym¬ 
metry was also used to determine Black Hole scattering in the extreme Kerr limit. As 
it happens, the correspondence seems exact, at least for a wide class of (vacuum and 
charged) Black Holes in four dimensions. 

While this seems surprising at first glance, one cannot help but wonder whether 
this fact is yet another facet on the intricate relationship between integrable systems 
and Conformal Field Theory [43]. General Relativity solutions like the (static, axisym- 
metric) black hole systems consist an integrable sector of Einstein equations in four 
dimensions [44], Requisition that the solutions of this integrable sector are singularity 
free basically imposes that there is a coordinate patch where the components of the 
metric are given in terms of rational functions. It follows then, apart from the non¬ 
trivial fact of separability of the equations, that the differential equations governing 
the propagation of fields in these backgrounds have its singularities easily classified, 
being Fuchsian or having a higher Poincare index. Moreover, the isomonodromy flow 
provides a non-linear symmetry of these equations stemming from the theory of flat 
holomorphic connections and the Painleve property. All in all, everything conspires to 
have the conformal blocks “solving” Fuchsian-type of differential equations, since con¬ 
formal blocks are in principle extracted from the representation theory of the Virasoro 
algebra. We can then describe the chain of deep mathematical facts pointing to the 
relationship between four-dimensional black holes and Conformal Field Theory. It is a 
very interesting open problem to see if this holds in the generic case. 
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